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Colliding Stellar Wind Models with Nonequilibrium 
Ionization: X-rays from WR 147 
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ABSTRACT 

The effects of nonequilibrium ionization are explicitly taken into account in a numerical 
model which describes colliding stellar winds (CSW) in massive binary sytems. This 
new model is used to analyze the most recent X-ray spectra of the WR+OB binary 
system WR 147 . The basic result is that it can adequately reproduce the observed 
X-ray emission (spectral shape, observed flux) but some adjustment in the stellar 
wind parameters is required. Namely, (i) the stellar wind velocities must be higher 
by a factor of 1.4 — 1.6; (ii) the mass loss must be reduced by a factor of ~ 2. The 
reduction factor for the mass loss is well within the uncertainties for this parameter 
in massive stars, but given the fact that the orbital parameters (e.g., inclination angle 
and eccentricity) are not well constrained for WR 147 , even smaller corrections to 
the mass loss might be sufficient. Only CSW models with nonequilibrium ionization 
and equal (or nearly equal) electron and ion postshock temperature are successful. 
Therefore, the analysis of the X-ray spectra of WR 147 provides evidence that the 
CSW shocks in this object must be collisionless. 

Key words: stars: individual: WR 147 - stars: Wolf-Rayet - X-rays: stars - shock 
waves 



1 INTRODUCTION 

The Wolf-Rayet (WR) star WR 147 (AS 431) is a clas- 
sical example of colliding stellar wind (CSW) binary: 
it possesses relatively strong X-ray emission and it is 
a nonthermal radio source. Probably the most interest- 
ing piece of information about this WR+OB binary sys- 
tem comes from the radio, namely, high-resolution inter- 
ferometer observations have spatially resolved its emis- 
sion into two components: a southern thermal source, 
WR 147S (the WN 8 component in the binary), and a 
northern nonthermal source, WR 147N (I Abbott et al.ll 19861: 
Moran et all 1 19891: Ichurchwell et all 1 19921: IContreras et all 



19961 : IWilliams et ail Il997i ; ISkinner et all 119991 . hereafter 
S99). Moreover, long-term radio variability has been re- 
porte d both from th e thermal |Contreras fc Rodriguez! 
1999) and nonthermal |Setia Gunawan et alj|200lh sources. 
The radio spectrum of WR 147S has a spectral index a = 
+0.6; +0.62 (S v oc v a ), which is very close to the canon- 
ical value for the thermal free-free e mission from ionized 
stellar winds (|Churchwell et ail 1 19921 : IWilliams et al. 1997; 
S99). On the other hand, the nonthermal radio emission of 
WR 147N cannot be described by a simple power law, which 
is expected from synchrotron models that assume power- 
law electron energy distribution (S99), and more elaborated 
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mode lling is required l|Doughertv et al J 120031 : 1 Pittard et ail 
l2006h . 

High-resolution infrared images clearly showed the NIR 
counterparts of the two radio sources located at nt 0'.'64 
from each other and WR 147N was classified as a B0.5V 
star (|Williams et al.|[l997l ) but an earlier spectral class, 08- 
09 V-III, wa s suggested from th e Hubble Space Telscope 
observations (|Niemela et ail Il998l). Given the d istance to 
WR 147 of 630 ± 70 pc jChurchwell et ailll992l ). the pro- 
jected (or minimum) binary separation is 403 ± 13 au. 
It is worth noting t hat analyzing radio and IR images, 
IWilliams et ail \ 19971 ) have found that the WR 147N radio 
peak is displaced by 0'.'07 southward from its IR counterpart. 
Although this value is within the uncertainties of the data, if 
it is real, then the nonthermal radio source can be associated 
with the CSW region located off the surface of the OB com- 
panion. This physical picture found further support from 
the high spatial resolution Chandra image, sh owing that the 
X-ray emission peak lies at or near WR 147N ijPittard et ail 

I200: 



The X-ray observations of WR 147 have revealed the 
presence of thermal emission from high temperature plasma 
but they also demonstrated how it is imporatnt to have data 
with good quality. Namely, WR 147 was first detected by the 
Eistein observatory and it was concluded tha t the X-rays are 
due t o plasma at temperature kT ^ 0.5 keV ijCaillault et ail 
Il985l ). Analysis of the moderate resolution CCD spectra, ob- 
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tained with ASCA and having not very high photon statis- 
tics, suggested that the X-rays likely originate in multitem- 
perature plasma whose cooler component has a temperature 
of kT fn 1 keV, but the characteristics of the hotter com- 
ponent could not be constrained (S99). Finally, the higher 
signal-to-noise XMM-Newton data of WR 147 presented by 
Skinner et al. (2007, hereafter S07) have undoubtedly re- 
vealed the elusive high temperature plasma in this object by 
detecting the Fe Ka complex at 6.67 keV. The correspond- 
ing analysis of the X-ray spectra has shown that the plasma 
tempeature is as high as kT = 2.7 keV. The authors pointed 
out that neither long-term X-ray variability nor changes in 
the X-ray absorption were detected. Thus, the 'evolution' 
of the plasma characteristics, mentioned above, is merely a 
result of the acquisition of higher quality data with each 
successive generation of X-ray satellites. 

Apart from the clear evidence for presence of high tem- 
perature plasma in WR 147 , S07 have found indications that 
further improvement of numerical CSW models is needed. 
Namely, higher terminal wind velocities are required to ac- 
count for the hot plasma detected in this object, and future 
CSW models must include some additional physics such as 
nonequilibrium ionization (NEI). 

Motivated by this, the aim of this study was to develop 
new CSW models that take into account the NEI effects and 
then to confront them with the high quality XMM-Newton 
spectra of WR 147 . The model details are given in § [3] The 
results from the model application are presented in § [4] and 
discussed in §[S] Our conclusions are found in §[15] 



2 DATA 

The X-ray observations of WR 147 were performed with the 
XMM-Newton observatory in November 2005. All the details 
about the EPIC data and the data reduction are found in 
S07. We only note that our analysis is based on the PN and 
MOS spectra of WR 147 and the latter is the sum of the 
spectra from the two MOS detectors. Thus, each analyzed 
spectrum has about the same number of counts (4436 in PN; 
4148 in combined MOS). They were rebinned to have a min- 
imum of 20 cts per bin. All spectral fits were performed in 
the recent version (11.3.2) of th e software pa ckage for anal- 
ysis of X-ray spectra, XSPEC |Arnaudlll996T l. The spectra 
were analyzed simultaneously with exactly the same model 
and only the XSPEC normalization parameter was decou- 
pled to allow for any calibration difference between the PN 
and MOS detectors. 



3 MODEL 

The goal of this study is to use a realistic model that deter- 
mines the physical parameters of the X-ray emitting plasma, 
and on its basis to simulate the X-ray spectra to be con- 
fronted with the observational data. A good starting point 
is to consider a simple model which describes the basic phys- 
ical situation and then elaborate on those that include more 
complex physics (e.g., nonuniform gas flows, thermal con- 
duction etc.) which might be potentially important for the 
case under consideration. As in the case of other wide stellar 



binaries (e.g., WR 140; IZhekov fe Skinner! [2000). an adia- 
batic model of CSW interaction was adopted for describing 
this phenomenon in WR 147 while more complicated models 
are beyond the scope of the present work. Adiabatic CSWs 
also have the advantage that the hydrodynamic solution in 
the interaction region does not depend on such details as: 
e.g., exact metal abundances, whether or not the plasma is 
in ionization equilibrium, or any difference between electron 
and ion temperatures. On the other hand, any of these may 
directly affect the shape of the emitted X-ray spectrum and 
the strength of various emission lines. Thus, it might be nec- 
essary to take them into account when simulating the X-ray 
emission from such objects. We note that all these physical 
characteristics allow the X-ray spectral simulations to be de- 
coupled from the hydrodynamic modelling of CSWs in wide 
binaries. Namely, the temperature and density distributions 
of the hot plasma, as initally derived from the numerical 
hydrodynamics, can be further postprocessed to model the 
corresponding X-ray emission. We next describe the steps in 
such an approach as applied to the case of WR 147 . 

3.1 CSW Model 

A primary step is to determine the distribution of tempera- 
ture and density of the X-ray emitting plasma. For this pur- 
pose, the hydrodynamic model of adiabatic CSWs by Myas- 
nikov & Zhekov (1993, hereafter MZh93) was used to derive 
the physical parameters in the interaction region of WR 147 . 
This model assumes spherical symmetry of the supersonic 
stellar winds that have attained their terminal velocities in 
front of the shocks. This assumption is well justified for wide 
stellar binaries where neither the radiative braking nor the 
orbital motion is expected to play an important role for the 
wind dynamics. In this case, the interaction region has cylin- 
drical symmetry and the gasdymanic problem is determined 
entirely by the stellar wind parameters and the binary sepa- 
ration. We note that the 'shock fitting' technique is used in 
this numerical model, and thanks to this an exact solution 
to all discontinuity surfaces (the two shocks and the con- 
tact discontinuity) is derived. This means that there is no 
numerical 'mixing' between the shocked gases of the stellar 
winds which further facilitates modelling the X-ray emission 
from the hot plasma, especially, when the different chemical 
composition is concerned. The procedure of how the effects 
of different electron an d ion temperaures can be taken into 
account is described in IZhekov fe Skinner! (|2000l ) and that 
for handling the nonequilibrium ionization is given below. 

3.2 NEI Effects 

In many astrophysical examples, the plasma emission is de- 
termined by the state of ionization equilibrium. If plasma 
characteristics such as temperature and density evolve on 
relatively small timescales, some time is also needed for a 
new state of ionization equilibrium to be established. There- 
fore, we might have a chance to witness transition phases 
in the plasma evolution. One such case is when the astro- 
physical plasmas are heated by shocks: gas temperature and 
density experience abrupt change. Fast shocks, as in super- 
nova remnants or CSW binaries, can result in high plasma 
temperatures, and the corresponding hot-plasma emission is 
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entirely determined by the collisional processes of ionization 
and excitation mostly with electrons. A detailed description 
of related processes that determine the hot-plasma emission 
both in and out of i onization equil ibriu m can be fou nd in 
the review papers bv lLiedahll l| 19991 ) and lMewel l|l999l ). 

Obviously, the transition phases, such as the state of 
nonequilibrium ionization, are important in evolving sys- 
tems. On the other hand, the NEI effects might play an 
important role also in steady state conditions as in the 
CSW shocks in wide binaries. Namely, despite the fact that 
the general picture does not change in time, each individ- 
ual parcel of gas follows the same 'evolutionary' pattern: 
it gets heated on the shock front, and then its ionization 
state relaxes downstream towards the corresponding equi- 
librium conditions: collisional ionization equilibrium (CIE). 
Whether or not the transition phase is important depends 
on how quickly CIE is established with respect to the typical 
timescale of the flow. 

In order to determine the relative importance of the NEI 
effects in CSW binaries, we consider a dimensionless param- 
eter T ne i = r/r ne i. Here r is the timescale of gasdynamics 
and T ne i is the characteristic time to establish ionization 
equilibrium, which we can refer to as the NEI timescale. 
The NEI effects must be taken into account if T ne i ^ 1 
but can be neglected for T ne i 2> 1. In the hydrodynamic 
model at hand, r = D/Voo, where 2D — a is the separa- 
tion between the binary components and Voo is the terminal 
stellar wind velocity. Unfortunately, there is no 'unique' NEI 
timescale and such a parameter can be introduced for each 
chemical element whose ionization state is considered. Al- 
though the latter requires one to solve a set of differential 
equations, qualitative considerations show that CIE is con- 
trolled mostly by the rate coefficients of collisional ionization 
and recombination to hydrogen- and helium-like species of a 
given chemical element. This is particulary the case for hot, 
X-ray emitting plasmas that are of interest in this study. 
Guided by this (see also the NEI discussion and the case 
example of oxygen in Licdahl 1999), we can introduce a rep- 
resentative NEI timescale r ne i = 10 12 /n e , where n e is the 
electron number density of the hot plasma. In that case T n€ .i 
can be expressed in terms of the basic CSW parameters: 



differential equations represents the evolution of the ioniza- 
tion state of each chemical element: 



—r— = Ci-iUi-i — (d + OLi)rii + Cti+iTli+i 



(2) 



Tnei — 1.21 



fiV? 000 D 16 



(1) 



where M5 is the mass-loss rate in units of 10 -5 M0 yr" 1 , 
V1000 is the stellar wind velocity in units of 1000 km s _1 , 
Die = -D/10 16 cm, fi, is the mean atomic weight for nucleons, 
and Xe = rie/n is the relative electron number density (n is 
the nucleon number density). 

Given the stellar wind and binary parameters of 
WR 147 (see § [31} , it is anticipated that the NEI dimen- 
sionless parameter F ne i « 1 and T ne i <C 1 for the shocked 
WR and O wind, respectively. This means that the NEI ef- 
fects likely play an important role for the X-ray emission 
from the CSW region in this object. 

Keeping in mind the cylindrical symmetry of the hydro- 
dynamic problem, a full set of 2D, time-dependent partial 
differential equations must be considered if the NEI effects 
need be taken into account in the CSW models. From a nu- 
mercial view point, an alternative to this is to consider the 
ionization balance along the flow streamlines. In this case 
the numerical task is much simpler and a set of ordinary 



where C and a denote the collisional ionization and recom- 
bination rate coefficients, respectively, and m-i,i t i+i are the 
number densities of any three consecutive ionization stages 
of a given chemical element. The parameter r e is the so- 
called ionization time, r e = J n e dt, and the integration 
starts at the shock front and is carried out downstream along 
a specific streamline. 

It is well known that the ionization equations (eqs. [2]) 
are a set of stiff ordinary differential equations, and vari- 
ous techniques have been developed to deal with this nu- 
merical problem. Very o ften an eigenfunction method (e.g. 
iHughes fc Helfandlll985T l is used to solve these equations effi- 
ciently, and this meth od was adopted in various NEI model s 
available in XSPEC l|Borkowski. Lverlv fc Reynolds! lioOlT l. 
Therefore, the latter was the basis for incorporating our 
NEI CSW models in the spectral fits to the X-ray data of 
WR 147 . We note that the NEI models in XSPEC usually 
assume constant electron density and temperature while this 
is not the case in the density and temperature stratified in- 
teraction region in CSW binaries. To overcome this prob- 
lem, we adopted a piecewise integration approach. Namely, 
the hydrodynamic model allows for deriving all the neces- 
sary plasma parameters (e.g., number density, electron tem- 
perature, ionization time) at the grid points along a given 
streamline of the flow. It is thus assumed that density and 
temperature are constant in between any two consecutive 
(j and j + 1) grid points along the streamline and now the 
NEI problem can be solved using the corresponding models 
available in XSPEC (with the only modification that at each 
step the initial values for eqs. [2] are those from the previous 
grid point). Once this is done, the X-ray spectrum of the hot 
plasma along a given streamline is derived and the sum over 
all the streamlines gives the total X-ray emission from the 
interaction region. 



3.3 CSW Models in XSPEC 

In order to make a direct comparison between the observa- 
tions and the X-ray emission predicted by the CSW models, 
we developed a new spectral model for XSPEC . The X- 
ray emission from the interaction region is thermal, thus, its 
characteristics are determined by the temperature, emission 
measure and chemical abundances of the hot plasma. This is 
true for the case of CIE and an additional parameter (ioniza- 
tion time) should be provided as well, if the NEI effects are 
important. All these necessary ingredients are given by the 
hydrodynamic CSW model. For modelling the X-ray emis- 
sion in CIE the XSPEC optically-thin plasma model apec is 
used at each grid point of the interaction region, while the 
case of NEI is considered as described above (§ 13. 2[) . Such 
an approach also has the advantage that the fits with the- 
oretical spectra of CSW binaries can be easily confronted 
with the results from the standard fits that usually make 
use of discrete-temperature plasma models (both in and out 
of ionization equilibrium) available in XSPEC (e.g., vmekal, 
vapec, vpshock). This, in turn, will be helpful for interpreting 
the X-ray data from CSW binaries when detailed numeri- 
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cal modelling is not feasible for various reasons (e.g., data 
quality, binary parameters are not available). 

Two types of test were performed for the new spec- 
tral models. First, for internal consistnecy, the distribution 
of emission measure of the hot plasma, as derived on the 
original grid of the hydrodynamic CSW model, was com- 
pared with the one resulting from the intergration along the 
flow streamlines. Second, results from our piecewise inte- 
gration approach were confronted with the ones from the 
XSPEC vpshock model. These tests made us confident about 
the further use of the spectral CSW models for analyzing 
the XMM-Newton data of WR 147 . 

Finally, using the known distance to the observed ob- 
ject, the emission measure from the hydrodynamic CSW 
model can be converted into XSPEC units 0. In such a case, 
the XSPEC norm parameter of our spectral model would 
have values close to unity if theoretical and observed fluxes 
match each other. Deviations from unity may be an indica- 
tion of inaccuraices in the adopted stellar wind parameters 
or the distance to the object. 



3.4 Basic Parameters for the CSW Model in 
WR 147 

For consistentcy with previous works (S99, S07), we define 
a 'standard' CSW model for WR 147 , which will also serve 
as a starting point for the fits to the X-ray emission from 
this object. As already mentioned, the distance of 630 pc 
to WR 147 is adopted in this study, and the stellar wind 
parameters are: M(WR) = 4 x 10" 5 M Q yr -1 , Voo{WR) = 
950 km s- 1 , M(0) = 6.6 x 10~ 6 M yr" 1 , v^O) = 
1600 km s" 1 . These values define a wind momentum ratio 
n = [M(0)v ao (0)}/[M(WR)v^(WR)] = 0.028 (note that 
A = 1/rj is used in MZh93 for the wind momentum ratio, 
but for consistency with the work of other authors r\ will be 
used throughout this study). Also, a minimum value of the 
stellar separation in this binary system is a — 2D — 403 au, 
and we note that this parameter depends on the inclination 
angle which is uncertain. Base d on analysis of VL A radio im - 
ages, values of i = 45° ± 15° dContreras fe Rodriguez! [l999l ) 
and i « 30° l|Contreras. Montes fe Wilkinl 120041 ) have been 
reported but other values may well be possible. It is so since 
the inclination angle was derived by fitting a thin shell model 
to the radio image and it should be kept in mind that such 
a model can not be considered satisfactory for wide stellar 
binaries, where the CSWs are adiabatic, since it gives only 
the shape of the contact discontinuity surface. On the other 
hand, the shape of the shock surfaces (e.g., their 'opening 
angle') is quite different which makes results from such a 
model application higly uncertain. Therefore, values for the 
binary separation in WR 147 a few times larger than the 
minimum value given above are not ruled out. 

Thus, we have all the necessary input parameters for the 
hydrodynamic CSW model at hand (MZh93), and elemental 
abundances have to be provided for modelling the thermal 



1 Conversion of emission measure into XSPEC units is done by 
multiplying it with the term 10~ 14 /Ami 2 , where d is the distance 
to the object in cm. This term comes from the XSPEC normal- 
ization factor of the flux for optically-thin plasma models. 



X-ray emission from the interaction region. Due to their ad- 
vanced evolutionary status, WR stars are objects with non- 
solar abundances and for this study the following set of them 
is adopted. The abundances for H, He, Ne and Ca are based 
on th e analysis of the IR spectrum of WR 147 l|Morris et alj 
2000), and t he other elements have values 'typi c al' fo r 
a WN star (|van der Hucht. Cassinelli fe Williams! Il986h . 
For convenience, all values are with respect to the solar 
number abundances (un ity means solar abundance as in 
lAnders fc Grevessdll989l ): H = 1, He = 25.6, C = 0.9, N 
= 140, O = 0.9, Ne = 44.3, Mg = 14.4, Si = 15.2, S = 3.6, 
Ar = 5, Ca = 4.7, Fe = 7, Ni = 1 (no information is available 
for Ni and it was assumed solar but we note that the spec- 
tral fits are not sensitive to its value). Due to the high X-ray 
absorption, the portion of the spectrum with good statistics 
(E ^ 0.9 keV) is sensitive only to the elemental abundances 
of Ne, Mg, Si, S, Ar, Ca, Fe, which are allowed to vary in 
our fits. In all fits the elemental abundances of the O wind 
were strictly solar and they are kept fixed in the modelling of 
the X-ray emission from WR 147 . The assumed WN abun- 
dances correspond to a relative electron density Xe = 1-71, 
mean particle weight (electrons and nucleons) /i = 1.16, and 
mean atomic weight for nucleons p, = 3.14. For the O star the 
adopted abundances give \e — 1-1, M = 0.62 and fi, — 1.3. 



4 RESULTS 

The results from previous application of the hydrody- 
namic CSW model for calculating the X-ray emission from 
WR 147 indicated that there is a luminosity mismatch be- 
tween the model predictions and the observations (S99, 
S07). Namely, the theoretical luminosity (or the amount, 
emission measure, of the hot plasma) is ~ 3 — 4 times 
larger than observed. Our new spectral analysis confirmed 
that this difference between theory and observations is 
now valid for both cases of CIE and NEI, if the nomi- 
nal wind parameters are used (e.g., the ones defined in 
§ 13. 4|) . This conclusion also holds if the mass loss of the 
O star, M(0), is varied. Two such cases were considered 
which have n = 0.02 (M(0) = 4.7 x 10~ 6 M yr" 1 ) 
and 0.04 (M(0) = 9.4 x 10~ 6 M Q yr" 1 ), respectively 
(note that all other parameters, M(WR), Voo (WR), Voo(0), 
kept their nominal values, § I3.4|l . Moreover, in all the three 
cases, r\ = 0.02, 0.028 and 0.04, the overall quality of the 
fits is not acceptable as indicated by the values of reduced 
xt = 1-63 — 1.93 (with degrees of freedom v = 332 for all 
the fits considered throughout this study). Therefore, it is 
conclusive that the CSW model not only predicts an amount 
of X-rays larger than that observed for WR 147 , but it also 
does not give a good match to the exact shape of the ob- 
served X-ray spectrum, especially, in its high-energy part 
above 4 keV, where the theoretical spectrum is softer than 
the observed one. Obviously, the data require presence of a 
hotter plasma in the interaction region and, as already sug- 
gested, a remedy to this problem could be the assumption of 
higher values for the stellar wind velocities (S07). Note that 
the contribution of the shocked O wind is no more than 12- 
15% at energies above 4 keV (S07 and this work), therefore, 
the velocity of the WN wind is a crucial parameter. 

To explore this possibility, a set of CSW models (both in 
CIE and NEI) was considered with increased wind velocities 
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Figure 1. Results from the CSW model fits. Left column: Mass-loss rates scaled while the binary separation is kept at its 'nominal' 
value. Right column: The binary separation is scaled while mass-loss rates have their 'standard' values. Upper panels: Reduced xt 
vs. scaling factor for the wind velocities. Dotted lines denote the models with nonequilibrium ionization (NEI), and dashed lines present 
models that assume the hot plasma is in collisonal ionization equilibrium (CIE). The plus signs are for the case of r\ = 0.04; asterisks 
mark models with rj = 0.028; and diamonds present models with rj = 0.02. The solid lines illustrate the effect of nonequal electron and ion 
temperatures Q3 = T e /T = 0.2) for r\ = 0.028 as triangles are for models with NEI while squares present those in CIE. Lower panels: 
The visual extinction as derived from the fits with various models where the Nh to Ay conversion formula from Gorcnstein (f975) was 
used. All symbols and lines have the same meaning. The two horizontal solid lines give the optical extinction deduced from IR studies 
(Ay = 11.5 mag, Churchwell et al. 1992; Ay = 11.2 mag, Morris et al. 2000). The degrees of freedom in all the fits are v = 332. 



for all the three basic cases (77 — 0.02, 0.028 and 0.04). It is 
anticipated that increasing wind velocities will make the NEI 
effects more pronounced (see eq. [I]). 

On the other hand, assumi ng the distance to WR 1 47 is 
well constrained (630 ± 70 pc, IChurchwell et al J 1 19921) the 
mass loss or the binary separation should be changed in 
order to solve the problem with the luminosity (emission 
measure) mismatch. Namely, it is necessary to decrease the 
emission measure in the interaction region which can be done 
in two ways: (i) by decreasing the mass loss rates; (ii) by in- 
creasing the binary separation. This follows from the relation 
between the emission measure in the interaction region and 
the parameters for spherically symmetric stellar winds: 

f M 2 

EM — I n^nudV oc — ~ — 

/ v ,,L) 

J wind 

where n e ,riH oc - — D z i an d V oc D 3 . As we see, increasing 
the wind velocity also means a smaller value for the emission 



measure, but we note that the general effect of this quan- 
tity is to control the shape of the resultant X-ray spectrum 
since the wind velocities are directly related to the plasma 
temperature in the interaction region. 

Thus, to discriminate between influences of various 
CSW parameters on the X-ray emission from WR 147 , the 
basic cases mentioned above, -q — 0.02, 0.028 and 0.04, were 
considered in two different ways. First, the mass loss rates 
of the both stellar winds were scaled simultaneously by a 
factor of 0.6 (77 = 0.02), 0.5 (77 = 0.028), and 0.4 ( V = 0.04), 
while the binary separation had its 'standard' value. Sec- 
ond, the mass losses kept their nominal values while the 
binary separation was increased correspondingly by a factor 
of 3 (tj = 0.02), 4 (77 = 0.028), and 6.25 (77 = 0.04). For each 
of the three basic cases, a grid of model fits was performed 
as the wind velocities changed with the same scaling factor 
that varied between 1 and 2. Adopting such an aprooach, 
namely, to scale the mass losses and/or the wind velocities 
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Figure 2. Examples of background-subtracted PN and MOS (MOS1+MOS2) spectra of WR 147 and the corresponding best fit NEI 
CSW models for r\ = 0.028 and stellar wind velocities scaled by a factor of 1.5. Left panel: the case of reduced mass losses by a factor 
of 0.5. Right panel: the case with increased binary separantion by a factor of 4. Fit results are given in Table [T] 



Table 1. Examples of CSW (rj = 0.028) spectral fits for WR 147 





Mass loss*"' 


Binary separation' 6 ' 


x7c!of 


379/332 


371/332 


N w (10 22 cm" 2 ) 


2.46 [2.30 - 2.64] 


2.48 [2.32 - 2.64] 


H 


1 


1 


He 


25.6 


25.6 


C 


0.9 


0.9 


N 


140.0 


140.0 





0.9 


0.9 


Ne 


16.9 [7.1 - 32.6] 


5.8 [0.3 - 13.7] 


Mg 


3.5 [1.9 - 5.9] 


1.6 [0.5 - 2.9] 


Si 


5.1 [4.1 - 6.3] 


3.7 [3.1 - 4.5] 


S 


7.8 [6.5 - 9.3] 


6.4 [5.4 - 7.4] 


Ar 


9.1 [5.4 - 12.9] 


8.5 [5.1 - 11.9] 


Ca 


9.4 [3.4 - 15.7] 


9.0 [2.8 - 15.3] 


Fe 


6.2 [4.6 - 7.9] 


6.9 [5.1 - 8.8] 


Ni 


1 


1 


F^' (PN) 


1.45 (20.1) 


1.45 (28.2) 


F£> (MOS) 


1.56 (21.6) 


1.56 (30.4) 



All abundances a re with respect to their solar values 
jAnders fe G rcvessc 1989). Brackets enclose 90% confidence in- 
tervals. 

a CSW model with decreased mass losses (0.5 of the nominal 
values), stellar wind velocities scaled by a factor of 1.5 (see text). 
6 CSW model with increased binary separation (by a factor of 4), 
stellar wind velocities scaled by a factor of 1.5 (see text). 
c The observed X-ray flux (0.5 - 10 keV) followed in parentheses 
by the unabsorbed value. Units are 10 — 12 ergs cm -2 s — 1 . 

of the both winds with the same factor has the advantage 
that the wind momentum ratio remains the same. This in 
turn facilitates the numerical hydrodynamic work needed 
prior to the spectral fits. 

Finally, it is worth noting that the decreased mass losses 
and the increased binary separation both have similar ef- 
fect, namely, they as well as the higher wind velocities make 
the NEI effects more pronounced (the F n ei values become 
smaller; see eq. [I]). 

Figure [1] presents the results from the fits to the X-ray 



spectrum of WR 147 for the cases discussed above. It is im- 
mediately seen that the CIE CSW models are not able to 
give acceptable fits (\t > 1-45). On the other hand, the in- 
clusion of NEI improves the quality of the spectral fits reach- 
ing acceptable values of xt = 1-10— 1.15, provided the wind 
velocities have values higher than the nominal ones by a fac- 
tor of 1.4—1.6. Interestingly, if the Gorenstein (1975) conver- 
sion formula is used, the derived X-ray absorption from the 
best spectral fits is in a good correspondence with the data 
for the visual extinction deduced from other studies. This 
fact is also encouraging for the conclusion that the develop- 
ments of the CSW model presented in this study are in the 
right direction. Examples of the most successful CSW mod- 
els are shown in Fig. [5] and details about the fits are given in 
Table [1] It is seen from the figure that the theoretical mod- 
els match very well the overall shape of the observed spectra 
with the exception of the soft energy part (E < 0.8 keV). 
We believe that future data with higher photon statistics 
will help resolve this issue. 

Thus, the results of the analysis of the X-ray emission 
from WR 147 show that there are at least three requirements 
for the numerical CSW model to be successful: (i) higher 
stellar wind velocities by a factor of 1.4 — 1.6; (ii) reduced 
mass losses for the stellar winds or larger binary separation 
(or both); (in) the NEI effects must be taken into account. 
We note that the needed reduction factor for the mass loss of 
~ 2 is well within the uncertainties for this parameter of the 
stellar winds in massive stars. But given the fact that the 
orbital parameters (e.g., inclination angle and eccentricity) 
are not well constrained for this binary system, much smaller 
corrections to the mass losses might be sufficient. 

On the other hand, the wind velocities in massive stars 
are constrained much tighter usually on the basis of the P 
Cygni profiles of ultraviolet and optical spe ctral lines. This is 
not ex actly the case of WR 147 . Namely, IChurchwell et al.l 
(1992) derived a wind velocity of 900 km s _1 from the 
analysis of the Hel 2.058pm P Cygni line profile, and the 
spectral resolution was 500 km s _1 . Also from analysis 
of NIR observations with spec tral resolution of 470 km 
s _1 (Hel 1.083/^m emission line). lEenens fe Williams! (|l994T ) 
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deduc ed 1100 km s for the wind veloc i ty in WR 147 . 
Later, lHamann. Koesterke fc Wessolowskil (|l995l) proposed 
a 1000 km s _1 wind velocity, as their analysis was based on 
optical Hel lines. Based on ISO-SWS data with spectral res- 
olution of 150-300 km s^ 1 , a wind velocity of 950 km s" 1 was 
derived from the line profiles of the fo rbidden lines [Nellll 
15.5/ j m, [SilV] 10.5Mm, [C alV] 3.21^m l|van der Hucht et all 
1 19961 ; iMorris et all [2000) . As we see, the wind velocity in- 
formation on WR 147 comes mostly from IR spectra which 
not always had a very high spectral resolution. For a large 
sample of galactic WR stars lEenens fc Williams! (| 1994! ) com- 
pared the terminal wind velocities derived from infrared data 
(Hel 1.083/im and 2.058/^m lines) with those from the con- 
ventional methods (based either on the P Cygni profiles of 
resonance lines in the ultraviolet region or on the widths 
of the optical lines). Interestingly, they have found that the 
values derived in IR are systematically lower by ~ 30%. 

If this is the case, then the wind velocity scaling factor 
derived here from the analysis of the X-ray data on the basis 
of the CSW model might not be unreasonable. Nevertheless, 
it is worth noting that if no solid evidence for higher stellar 
wind velocities are found in the future, this will pose serious 
problems for the CSW scenario as a central mechanism for 
the origin of X-rays in WR 147 . 



5 DISCUSSION 

5.1 CSW Shocks: Collisional or Collisionless? 

All the cases considered so far in this study assume a single- 
temperature hydrodynamic approximation, that is equal 
electron and ion temperatures. We recall that the thermal 
X-ray emission traces only the hot electrons, so, its anal- 
ysis gives information mostly on the electron temperature 
in hot plasmas. If the plasma is heated through shocks, the 
kinetic energy of the gas flow is immediately (at the shock 
front) deposited to heavy particles (nucleons) and energy ex- 
change between heavy and light particles (electron heating) 
occurs downstream via Coulomb collisions. Thus, some time 
is needed for the electron and ion temperatures to equili- 
brate. This timescale depends basically on the electron tem- 
perature and the nucleon number density ( Braginskii 1 19581 ; 
Spitzcr 1962), as it is longer in plasmas with higher temper- 
ature and lower density. 

Although the electron heating mechanism at the shock 
front is not yet fully understood, it should be kept in mind 
that there are two general cases: electron heating in colli- 
sional and collisionless shocks. In the former case, the elec- 
trons initially attain a very low temperature mainly through 
adiabatic heating (Zeldovich & Raizer 1967). The situation 
could be very different in the latter case when various heat- 
ing mechanisms might be at work and a range of degrees 
of temperature equilibration is possible up to equal elec- 
tron and ion temperatures (e.g..lMcKeelll974l ; |Papadopoulus! 
Il988l ; ICargill fc Papadopouluslll988l '). These possibilities are 
usually described by a simple parameter that relates the 
electron temperature at the shock front to the postshock 
mean plasma temperature, f3 = T e /T (thus f) m 0. — 1.). 

As the effects of unequal electron and ion tempera- 
tures are important in high temperature, rarefied plasmas, 
they are thus expected to play a role for colliding stellar 



wi nds in wide binaries. Suc h models were first considered 
bv lZhekov fc Skinner! l|2000h who also introduced a dimen- 
sionless parameter that determines the relative importance 
of electron heating via Coulomb collisions downstream in the 
interaction region (see eq. [1] in their work). For the nominal 
wind and binary parameters of WR 147 (§ 13. 4p . the value of 
this dimensioneless quantity suggests that the temperature 
equilibration effects can be neglected when modelling the 
X-rays from this object. 

On the other hand, as was shown in § [4] the success- 
ful CSW models require higher wind velocities, and reduced 
mass loss or larger binary separation than the nominal val- 
ues of these parameters. This in turn means a slower tem- 
perature equilization, that is a higher chance for the unequal 
electron and ion temperatures to affect the X-ray emission 
from WR 147 . To explore this possibility, we ran models 
(both NEI and CIE) with (3 = 0.2, and the main result 
is the lower quality of the spectral fits (see Fig. [T|. Obvi- 
ously, the reson for this is that the temperature equalization 
timescale is relatively long, therefore, the electron temper- 
ature remains below the mean plasma tempearture in the 
interaction region. Guided by this and the fact that models 
with CIE plasma cannot reproduce the X-ray spectrum well, 
more CSW models with different values of were consid- 
ered as preference was given to cases with velocity scaling 
factor larger than 1.5. Some results from the spectral fits 
with these NEI models are shown in Fig. [3] It is seen that 
acceptable fits (xt = 1-10 — 1.15) are possible if the wind 
velocities are of 1.5 — 2.0 times larger than their nominal 
values, and for a high degree of the temperature equilization 
(P 0.4). The best fits for high velocity scaling factors are 
obtained with a lower initial electron temperature relative 
to the mean postshock gas temperature. 

It is then interesting to check whether good spectral fits 
are possible also for very low values of the initial electron 
temperature. To do so we ran NEI CSW models with f3 = 
0.001, a value that could be considered typical for collisional 
shocks. The results are shown in Fig. U and it is seen that 
no acceptable fits are possible for the 'reasonable' values of 
the velocity scaling factor. 

Therefore, it is conclusive that the CSW shocks in 
WR 147 are collisionlessl 



5.2 CSW Models vs. Discrete Temperature 
Models 

The analysis presented here has shown that the CSW model 
can explain the XMM-Newtion X-ray data for WR 147 , 
provided that some reasonable adjustments to the currently 
accepted mass-loss parameters and orbital separation are 
made. On the other hand, it was shown in S07 that much 
simpler models, discrete temperature ones, could be equally 
successful, as the quality of the spectral fits shows. The obvi- 
ous advantage of the CSW model is that it describes a real 
physical situation which is likely present in massive bina- 
ries, but what could we learn from the application of these 
simpler models? 

It is worth recalling that the interaction region in CSW 
binaries is temperature stratified, thus, its X-ray emission 
is a sum of individual contributions from plasmas with a 
range of temperatures weighted by their emission measure. 
It is then not surprising that in such a case the total ('in- 
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Figure 3. Quality of the fit for NEI CSW models (r] = 0.028) with different degrees of initial temperature equilibration (J3 = T e /T). 
Shown are cases with wind velocity increased by a factor of 1.5, 1.7, and 2.0, respectively. Left panel: Mass-loss rates scaled while the 
binary separation is kept at its 'nominal' value. Right panel: The binary separation is scaled while mass-loss rates have their 'standard' 
values. The degrees of freedom in all the fits are v = 332. 



tegrated') X-ray spectrum can be quite well represented by 
isothermal plasma emission with some 'average' tempera- 
ture, and its value must be lower than the maximum plasma 
temperature in the shocked plasma. We note that the suc- 
cessful NEI CSW models require a wind velocity scaling fac- 
tor of at least ~ 1.4 (see § . Given the nominal wind pa- 
rameters (§ I3.4|l . this suggests for the shocked WR wind a 
maximum plasma temperature of ~ 4 keV, which is found 
near the line of centers in CSW binaries (for a strong shock 
with adiabatic index 7 = 5/3 the postshock temperature is 
kT = 1.956/xdiqqq , keV, where viooo is the shock velocity 
in units of 1000 km s _1 ). From such a point of view it is 
understandable why a constant temperature plane-parallel 
shock model with kT ~ 2.7 keV does a good job (S07). This 
model (vpshock in XSPEC ) has the most important char- 
acteristic of the CSW models for WR 147 , namely, it takes 
into account the NEI effects. As also discussed in S07, the 
fact that the temperature of this shock model was higher 
than the maximum expected for the currently accepted WR 
wind velocity (950 km s _1 ) was already indicative that 
higher stellar wind velocities might be present in this ob- 
ject. These authors reached similar conclusions on the basis 
of their discrete two-temperature models with CIE. 

Thus, the use of simpler models for representing the X- 
ray emission from hot plasmas can, on the one hand, be very 
helpful by guiding further modelling that elaborates a real 
physical picture in the studied objects. On the other hand, 
it may also pose some problems for interpreting the results 
when models bearing quite different physics (e.g. CIE vs. 
NEI) are equally successful, as is the case for WR 147 . A 
solid argument in favour of the NEI case here could be that 
the more elaborated models are successful only if these ef- 
fects are taken into account. But, it will be much more valu- 
able if a purely observational evidence is found that helps 
resolve this issue. For example, future gratings data with 
good photon statistics will detect relatively strong forbid- 
den lines in the He-like triplets of various elements if the 
NEI effects play an important role. Such observational data 
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will thus reveal where exactly the X-ray emission forms, i.e. 
in the rarified hot plasma of the CSW interaction region. 

5.3 CSW Models and Nonthermal Emission from 
WR 147 

We recall that wide CSW binaries are not only sources of 
strong X-ray emission but are non-thermal radio emitters 
as well, and WR 147 is no exception to the rule. Rela- 
tivistic electrons must be available in such objects, thus, 
the strong shocks that shape the CSW zone in these sys- 
tems a re likely places wher e the electrons are being accel- 
erated. iPittard et al.l l|2006l) have presented this idea in a 
quantitative manner developing a numerical model of the 
synchrotron radio emission from CSW binaries. Although 
this model is not self-consistent it was the first attempt to 
consider the non-thermal radio emission from such objects 
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in detail. Naturally, their theoretical predictions were then 
confronted with the detailed radio observations of WR 147 . 
One of the main conclusions from this analysis was that 
the acceleration of relativistic electrons appear to be very 
efficient requiring up to half of the thermal energy of the 
postshock plasma to be deposited into non-thermal parti- 
cles. We must note that if this were the case, one of the 
basic consequences for the CSW zone would be a noticeable 
decrease of the temperature of the shocked hot gas. 

Contrary to this, our analysis of the XMM-Newton data 
for WR 147 has shown that the X-ray emitting plasma is 
hotter than suggested by the nominal parameters of the stel- 
lar winds in this object. It is then conclusive that there is 
evidence that the CSW model is capable of explaining the 
X-ray emission in the wide WR+OB binary WR 147 but 
further refinements are needed to successfully describe the 
corresponding non-thermal radio emission from this object. 



reveal strong forbidden lines in the He-like triplets of various 
elements. The detection of strong / lines would be the most 
solid argument that the hot plasma in WR 147 is in a state 
of nonequilibrium ionization. It will then be evident that the 
X-ray emission in WR 147 originates in rarefied plasma that 
is characteristic of CSWs in wide binary systems. 
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6 CONCLUSIONS 

We developed an improved colliding stellar wind model that 
includes the effects of nonequilibrium ionization in hot plas- 
mas, and the XMM-Newton X-ray spectra of WR 147 were 
analyzed in the framework of this model. The basic conclu- 
sions are the following. 

1. If the 'standard' (§ I3.4|l stellar wind parameters and 
binary separation are assumed, the CSW model predicts 
an emission measure of the hot plasma that is 3 — 4 times 
larger than what is needed to explain X-tay emission from 
WR 147 . Also, the shape of the theoretical spectrum does 
not match well that of the observed spectra, especially, at 
higher energies (E ^ 4 keV). It is then concluded that hot- 
ter plasma must be present in the CSW zone of this object, 
and as proposed by S07, such a hotter plasma could be the 
result of higher wind velocities. 

2. By exploring a range of stellar wind velocities, it 
is shown that CSW models with nonequilibrium ionization 
are able to satisfactorily fit the observed X-ray spectra of 
WR 147 , provided the stellar winds are faster by a factor 
of 1.4 — 1.6. To resolve the issue of the increased emission 
measure, the mass loss rate should be decreased by a factor 
~ 2 or the binary separation must be increased by a factor 
~ 4 compared to the nominal values of these parameters. All 
these parameter changes are within their corresponding un- 
certainties. It is worth emphasizing that in the framework of 
the CSW picture only NEI models can correctly represent 
the observed X-ray data from this object. 

3. Different electron and ion temperatures are usually 
found in shock heated rarefied plasmas, which is likely the 
case in wide CSW binaries. Models that assume various de- 
grees of temperature equalization in the CSW shocks were 
confronted with the X-ray data of WR 147 . Only NEI mod- 
els with very high temperature equalization were successful, 
a case typical for collisionless shocks. On the other hand, 
models that assume a very low electron postshock temper- 
ature, which is characteristic of collisional shocks, cannot 
explain the observed data. Therefore, the analysis of the X- 
ray spectra of WR 147 finds evidence that the CSW shocks 
in this object must be collisionless. 

4. If future X-ray grating observations of WR 147 with 
good photon statistics are obtained, they are expected to 
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